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Approximately universal optimality over several 
dynamic and non-dynamic cooperative diversity 
schemes for wireless networks 

Petros Elia and P. Vijay Kumar 


Abstract —In this work we explicitly provide the first ever 
optimal, with respect to the Zheng-Tse diversity multiplexing 
gain (D-MG) tradeoff, cooperative diversity schemes for wireless 
relay networks. The schemes are based on variants of perfect 
space-time codes and are optimal for any number of users and 
all statistically symmetric (and in some cases, asymmetric) fading 
distributions. 

We deduce that, with respect to the D-MG tradeoff, channel 
knowledge at the intermediate relays and infinite delay are 
unnecessary. We also show that the non-dynamic selection decode 
and forward strategy, the non-dynamic amplify and forward, 
the non-dynamic receive and forward, the dynamic amplify 
and forward and the dynamic receive and forward cooperative 
diversity strategies allow for exactly the same D-MG optimal 
performance. 

I. Introduction 

The emerging need for reliable communications of large 
quantities of data, at high rates, between small and inde¬ 
pendent users with small power supplies, no antenna arrays 
and minimal computational capabilities, brought to the fore 
cooperative-diversity in wireless networks, where distributed 
users relay messages for one another, in order to combat the 
fading and additive noise that hinters their joint communica¬ 
tion, hence improving the overall quality of service. Cooper¬ 
ation can be achieved by relating each user with a segment 
of a point-to-point communication scheme and essentially 
having the intermediate relays manipulate the signal from the 
information source in a way that the received signal at the final 
destination relates to that of a point to point transmission with 
multiple transmit antennas. Network outage analysis provides 
for the fundamental limits of the network’s performance. 

A. Structure of the paper 

In Section 2]we will describe the cooperative-diversity setup 
and the related notation that will be used in the rest of the 
paper, will briefly introduce the existing cooperative diversity 
strategies, state the existing performance bounds and quickly 
go over the basic aspects of the D-MG tradeoff. 

In Section [H] Theorem^ we will present the expression for 
optimal D-MG performance for the non-dynamic selection- 
decode-and-forward strategy (ND-SDAF [3], [4]) and will 
provide the exact coding methodology that achieves this 
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optimality. In proving the scheme’s optimality, we will use the 
existence of sets of approximately universal [7] codes which 
hold common elements, and which maintain approximate 
universality even when their structure is altered. We will 
also generalize with respect to network topology and fading 
statistics. 

In Section [TO] we will present a variant of the non-dynamic 
linear-processing relay network [1], Unlike with ND-SDAF, 
this non-dynamic receive-and-forward (ND-RAF) scheme re¬ 
quires for coding distribution across the relays. We will 
then proceed to analyze the second stage equivalent (‘two- 
product’) channel and prove its partial D-MG equivalence to 
the Rayleigh fading channel. Using the existence of sets of 
approximately universal codes with joint elements, we will 
present in Theorem |3 the optimal performance for the ND- 
AAF and ND-RAF scheme and the exact coding methodology 
that achieves this performance. As a guideline for other poten¬ 
tial constructions, we will introduce one by one all the neces¬ 
sary conditions for achieving optimality. A practical scheme, 
more suited for a network with a large number of users and 
which provides for a plethora of practical advantages, will 
also presented. Finally, a base station network setup will be 
presented together with some new variants of perfect codes 
that have the potential to render cooperation fruitful at the 
lowest possible SNR. 

In Theorem 0 of Section EY] we present D-MG optimal 
schemes for the dynamic amplify-and-forward and the dy¬ 
namic receive-and-forward strategies (D-AAF, D-RAF [5]). 
Optimality holds for any number of users, any set of channel 
statistics and any statistical asymmetry. 

Section offers a comparison of the cooperative diversity 
schemes, and through Theorem 0] it is concluded that all 
the above strategies have high-SNR outage regions of the 
same volume and thus achieve the same D-MG optimal 
performance. We point out that given statistical symmetry, 
the ND-RAF scheme offers this same optimal performance 
at a reduced delay, reduced decoding complexity, reduced 
signalling complexity, minimal computational efficiency and 
the highest ease for network deployment. Section IvTI presents 
some simulations. The rest of the sections are appendices. 

We begin we a general description of the network. 

B. Describing the network 

In [3], the authors describe the case where a set 

1Z = {Ri, R.2, ■ ■ ■ , Rm Rn+l) 


of n + 1 different terminals/relays, cooperate in their effort 
to communicate with each other. Each relay has the ability 
to communicate over n + 1 different orthogonal frequencies 
T = {u 1 , 1 * 2 ,-•' ,Vn,v n +i}. A certain relay Ri, wanting 
to communicate with relay d(Ri), broadcasts its original 
information over frequency v,. Depending on the availability 
of each intermediate relay, the set 

D(i? i )c{7e\{ J R i Ud(i? i )}} (1) 

is then the set of all intermediate relays that cooperate with 
Ri. Consequently, each relay Rj £ I) ( II,) transmits a possibly 
modified version of the received signal over frequency Vj. By 
the end of the transmission, d(Ri) has received the information 
from Ri over frequency u,, in a form of a superposition of 
faded versions of signals originating from Ri and from D(Ri). 

We will focus on the case where communication takes place 
in the presence of additive receiver noise, and in the presence 
of spatially independent quasi-static fading. Furthermore, we 
will assume complete knowledge of the fading channel at the 
receiver of the final destination, and depending on the cooper¬ 
ative diversity strategy, we will assume complete knowledge 
or absolutely no knowledge of the channel at the receivers 
of the intermediate relays. Finally, the half-duplex condition 
is imposed, due to practical considerations such as the large 
ratio between the transmission and reception powers at the 
relay antennas ([3], [4], [5]). 

a) Instance of a network: From [3] we see that without 
loss of generality we can analyze the overall network perfor¬ 
mance just be focusing on a snapshot of the network, as shown 
in Figure [T| where S is now the information source, I) the 
destination, f? 2 , • ■ • , R n are the intermediate relays, g, is the 



Fig. 1. Snapshot of wireless network, where terminal S utilizes its peers 
(i? 2 , • • • , Rn) for communicating with D. 

fading coefficient between S and intermediate relay II,. hi is 
the fading from II, to D and h\ is the fading coefficient from 
S to D. We consider h,, g, to be independently distributed 
circularly symmetric CA/”(0,1) random variables, remaining 
constant throughout the transmission. Vectors v t and w contain 
the elements Vi j and Wj corresponding to the additive receiver 
noise respectively affecting Ri and D at time t = j. Unless we 
state otherwise, we ask that all Vij and Wj be independently 
distributed CA/”(0,1) random variables. SNR will represent the 
ratio of the signal power to the variance of the noise at the 
receiver of D. 


C. Existing cooperative diversity strategies 

We proceed to briefly introduce the above mentioned coop¬ 
erative diversity strategies. 

1) Non-dynamic selection-decode-and-forward: In [4], 
Faneman, Tse and Wornell, described the non-dynamic 
selection-decode-and-forward cooperation strategy for the 
three-node network which asks for f ?2 to cooperate, by first 
decoding and then re-encoding, if and only if it is not in 
outage with respect to the information source S. The strategy 
requires full channel knowledge at the receivers of relay If 
and destination D. Focusing on the case where each node has 
a single transmit-receive antenna and where the fading scalars 
are from the Rayleigh distribution, the scheme has f? 2 decode 
if and only if 

Source <l0g 2 (l + SNR| ff2 | 2 ) (2) 

where R somC e is the information rate of the transmission of S, 
measured in bits per channel use (bpcu). Again specific to the 
case of Rayleigh fading and of having a single transmit-receive 
antenna per node operating in the half-duplex environment, it 
was shown in [4] that given a source-to-destination informa¬ 
tion rate R , measured in bits per network channel use (bpncu), 
and given a multiplexing gain [6] 

R 

r ~ log 2 (SNR)’ 

then the high-SNR probability of network outage, and thus the 
optimal D-MG performance, is given by: 

j ( ,\ -^outage (f) _ on „ x _ 1 /o 

^ND-SDAFoutv ) •— i CXTD — 2(1 2r), r m ax — 1/2 

log SNR 

and that this optimal performance will be achieved by some 
random Gaussian codes of infinite length. This result is readily 
generalized to a network of n + 1 users where the cooperation 
strategy for intermediate relay Rj is dictated by 

Rj£D(S) f?<log 2 (l+SNR| 5j | 2 ), j = 2, •■•,». (3) 

Furthermore, given the implicit knowledge of the multiplexing 
gain at the nodes 1 , and conditioned on infinite time duration, 
the same approach translates to an optimal performance of 

rfND-SDAF-opt(r) = (n - 1) ^1 - + (1 - r) (4) 

where a + := max(0, a). The same performance was predicted 
in [4] for the non-dynamic amplify and forward scheme. 

Fetting r coop describe the multiplexing gain corresponding 
to the set of (f?,SNR) pairs where the D-MG performance 
of the cooperative scheme stops being equal to the D-MG 
performance in the non-cooperative case, we note that given 
a rate R, cooperation essentially applies only for SNR values 
greater than 

SNR > 2 r coop — 2 n — 1 R 

As in [6], we have =, > and < denoting asymptotic expo¬ 
nential equality and inequalities respectively. 

1 Knowledge of the code/rate and of the power constraint, together with the 
existing assumption of unit variance additive noise, jointly imply knowledge 
of the multiplexing gain. 
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Fig. 2. Existing D-MG bound for the single intermediate relay, ND-SDAF 
strategy. Given a network rate R (bpncu), cooperation applies for SNR > 



Fig. 3. The PEP based D-MG bound for the single intermediate relay ND- 
RAF strategy, predicts that given finite duration random coding, cooperation 
applies for SNR > 2 kR , k >> 1 and that given optimal coding, cooperation 
applies for SNR > 2 7R . 


2) Non-dynamic receive-and-forward : In a variant of the 
two-stage wireless relay network model proposed by Jing 
and Hassibi in [1], [2], we will consider the case where the 
n—1 intermediate relays are only allowed to perform linear- 
processing (time-averaging based on space-time codes) on the 
received signal. Knowledge of the channel ( hi,gi ) is given 
only to the receiver of the final destination D. 

In [1], [2], the authors also present a bound on the net¬ 
work’s pairwise error probability (PEP), which in the high- 
SNR regime is a function only of SNR and of the minimum 
eigenvalue of any difference of any two code matrices of 
the corresponding distributed space-time code which performs 
the linear-processing (see Appendix II1-0 . This result, in 
conjunction with the finite duration random coding proposed 
in [1], guarantees only for full-diversity but does not guarantee 
any bound on the code eigenvalues. As we discuss in Appendix 
IH-CI this implies that with finite delay random coding, the 
maximum achieved multiplexing gain can get arbitrarily small 

^nd-raf(t) > (n - 1)(1 - kr) + + (1 - r), k » 1 

thus potentially requiring for asymptotically high SNR, in 
order for cooperation to apply. 

A better lower bound on the optimal D-MG performance of 
the network is obtained when the same PEP result of [1], is 
applied on codes that are approximately universal and thus D- 
MG optimal over any channel. In the same Appendix III-CI we 
show that given such codes, expanding the existing PEP bound 
towards D-MG, guarantees a lower-bound on the network’s D- 
MG performance of: 

/ 4n-l\ + 

6^nd-raf( f) > (n - 1) f 1 - rc _ 1 r j + (1 - r). (5) 

This implies a maximum diversity d(0) = n and cooperation 
after 

SNR min > 2^ r . (6) 

3) Dynamic amplify-and-forward and dynamic receive-and- 
forward: In this D-AAF strategy, first introduced in [5], the 
original source S transmits at each time instance and each 
intermediate relay takes turns in transmitting an amplified 


version of a previously received signal, and do so only at even 
time indexes. Again in [5], it is concluded that for n+1 users, 
given infinite time duration, the optimal tradeoff is given as 

<A>aaf(V) = (n ~ 1)(1 - 2r) + + (1 - r). (7) 

From the analysis in [5], we can see that the same optimal 
D-MG result holds when the intermediate relays do not use 
channel knowledge. 

Having briefly introduced the cooperative strategies, we 
conclude the introduction with a quick exposition of the basic 
tools to be used in the network analysis and encoding. 

I). Preliminaries on the D-MG Tradeoff and approximate 
universality 

Let an n x T space-time code X operate at rate 
R x = y log 2 (|A’|) bpcu, and let r x be the multiplexing gain 
(normalized rate) given by 

R x 

r * = l°g 2 (SNR) 

corresponding to 

\X\ = SNR r * T . (8) 

For large SNR, the capacity over an n x n r Rayleigh fad¬ 
ing channel is given by C x ~ min{n, n r } log 2 (SNR), 
implying a maximum achievable multiplexing gain of 
max = min{n, n r }. The diversity gain corresponding to a 
given r x , is defined by 

us i • log(-Pe) 

d{r * } = -*te»bi(SN5)' 

where P e denotes the probability of codeword error. In a recent 
landmark paper, Zheng and Tse [6] showed that there exists a 
fundamental tradeoff between diversity and multiplexing gain, 
referred to as the diversity-multiplexing gain (D-MG) tradeoff. 
For a fixed integer multiplexing gain r x , and T > n + n r — 1, 
the maximum achievable diversity gain d(r x ) is shown to be 

d{r x ) = (n — r x )(n r — r x ). 


(9) 













The function for non-integral values is obtained through 
straight-line interpolation. For T < n + n r — 1 only bounds 
on the maximum possible d(r x ) were available [6]. It was 
also shown in [6] that there exist random Gaussian codes with 
T > n + n r — 1 that achieve the above optimal tradeoff. For 
such optimal codes, the probability of error coincides with the 
probability of outage 

P(outage) = P (error-optimal) 

The Zheng-Tse result sparked considerable interest in meet¬ 
ing this new D-MG frontier and explicitly providing D-MG 
optimal schemes. 

The problem of constructing explicit D-MG optimal ST 
codes over the Rayleigh-fading channel for any pair (n, n r ) 
was settled in [9] where it was shown that cyclic-division- 
algebra-based space-time codes having a certain non-vanishing 
determinant (NVD) property are optimal with respect to the 
D-MG tradeoff of the Rayleigh-fading channel. The same 
authors used this result to establish the D-MG optimality of 
constructions of space-time codes found in Belfiore et. al. 
[20] and Kiran and Rajan [22]. These prior constructions 
were restrictive in terms of the values of n that can be 
accommodated. In [9], [11], a general construction of D-MG 
optimal codes was provided that was valid for all n, n r and 
all T > n. 

Furthermore, as is shown in [7] (see also [8]), the above 
codes are approximately universal [7] and thus D-MG optimal 
over all slow-fading channels, independent of the channels’ 
statistics. As shown in the proof of Theorem 4.1 of [7] 
(Appendix A.2), due to the bound on their smallest eigenvalue, 
such codes satisfy the extra property that 

P(error | no outage) = e~ SNR+ = SNR~°° = 0. (10) 

In a slight abuse of notation, we will use the term ‘approx¬ 
imately universal given statistical symmetry’ to denote a D- 
MG optimal scheme that also satisfies 03, and does so for 
all channel statistics with identically distributed path fading. 

This concludes the introduction and we can now proceed 
with analysis and encoding for the first cooperative strategy. 

II. Optimality in the Non-Dynamic 
Selection-Decode-and-Forward Strategy 

In regards to the ND-SDAF strategy, we will here improve 
the existing optimal performance bound in 0 , generalize to 
a new optimality expression that holds for a larger family of 
channel statistics and explicitly construct the first ever coding 
scheme that guarantees for optimality in finite and minimum 
time delay, finite decoding complexity and for any number 
of users. We finally bound the performance for a variety of 
network topologies and explicitly provide encoding schemes 
that guarantee these bounds for any set of statistics. 

a) Cooperation strategy and general coding require¬ 
ments: We begin by noting that in this scheme, cooperation 
does not provide gains unless the probability of erroneous 
decoding at the intermediate relays is much smaller than 
the probability that some relay does not select to cooperate. 
Necessary for D-MG optimality is a decoding strategy that 


relates to outage as in 0 , in which case decoding at the node 
(during the first stage) is required to provide for 

P (error at Ri j no outage) << P (error at D with cooperation) 

which in terms of approximate universality, translates to 

P(error at Ri | no outage) = SNR~°° = 0. (11) 

D-MG optimal codes that guarantee that 

P(error) = P(outage) (12) 

can be found in the family of random Gaussian codes, which 
also maximize the mutual information of both the first and 
second stages of the transmission, and hence optimize the 
overall D-MG performance of an infinite duration network. 
The lack of eigenvalue bounds in the finite length version of 
these random codes does not allow for CD to hold, and instead 
as shown in [6], finite length random Gaussian codes satisfy, 

P(error | no outage) = P(outage) > SNR -00 . 

Consequently, the use of random Gaussian codes in the first 
stage of the network transmission (source-to-relay) immedi¬ 
ately imposes a requirement for infinite time duration [3], [4], 

[5]. 

We will proceed to explicitly provide a network encoding 
scheme that meets the outage region of any statistically sym¬ 
metric ND-SDAF network and does so in finite time duration. 

A. Explicit D-MG optimal encoding for the statistically sym¬ 
metric ND-SDAF network 

To achieve optimality, we jointly treat the first and second 
stage encoding methods by first constructing alxn approxi¬ 
mately universal code over SISO channels and then annxn 
approximately universal code over MISO channels, which has 
the extra two properties that it maintains its optimality even if 
it is truncated and that it has the optimal number of common 
entries with the first stage SISO code. 

b) First stage transmission and the horizontally- 
restricted, approximately universal, perfect code: For a relay- 
network where all the nodes have one transmit/receive antenna 
operating in half-duplex over a statistically symmetric channel, 
the proposed coding scheme asks for the source S to sequen¬ 
tially transmit, during time t = 1,2,--- , n, the n-length vector 
k = 9z where z comes from the 1 x n horizontally-restricted 
perfect code Xh, 

X h = {z=[z W * (2) z (n) ]} (13) 

= {[/l h ■ ' • fn]M n , V[/l h ■■■ fn] G A n } 

where the /, are from a discrete information constellation A 
such as QAM or HEX, that scales with SNR as \A\ = SNR r , 
and where M n is the unitary lattice generator matrix for perfect 
codes [21], [8], [10]. Finally 6 is the normalization factor 
such that E[|6G^ | 2 ]=SNR. The approximate universality of 
X h over the first stage SISO channel and the subsequent 
satisfaction of CD, are established by observing that the code 
carries the same information and has the same non-vanishing 
product distance (Xh consists entirely of one of the layers of 


the CDA-perfect codes) as the simple QAM scheme, which 
was shown in [7, Section 3] to be approximately universal 
over all SISO channels. 

c) Second stage transmission and the residual approxi¬ 
mate universality of CDA codes: By t = n, due to the approx¬ 
imate universality of X^, each intermediate relay If £ D(S) 
has correctly decoded k and will participate in the second 
stage of cooperation which will take place if and only if 
r < i (this choice of r will become clearer later on), in 
which case the network encoding scheme asks from each 
Ri £ D(S ), i = 2, • • • ,n to transmit feW (the i th element of 
k) at time t = n + i — 1, thus allowing the decoder of D to 
‘see’ a possibly truncated version of the diagonal restricted 
perfect code, Xd, given as: 


X d = {diag(z) = diag(f_ ■ M n ), V/ £ A n }. (14) 


As with X h , the corresponding information alphabet A is 
discrete and M n , as in d 1 3k represents the orthogonal matrix 
that generates a lattice in the maximal field L of the division 
algebra of the perfect codes (see equation J59H . The above 
choices make Xd approximately universal over any i.i.d. 
MISO channel, since the product-distance of any difference 
of diagonals in the matrices of Xd, 

n— 1 n— 1 / n ~ 1 \ /n—1 \ 

n 5 z w =n ^ (n fkPk ) — -^L/Z[i] f ^ ^ fkPk j 

j =0 j =0 ^k =0 ' ^k =0 ' 

is an algebraic norm and thus a Gaussian integer, meaning that 


n—l 


n 6z d) 

3=0 


> 1 , 


which together with the fact that \A\ = SNR’’, satisfy all the 
related conditions in [7, Theorem 4.1], 

The above approximate universality only relates to the event 
where all intermediate relays Ri are in D(S). As we have 
seen though, it is the case that some relays might be unable 
to decode, in which case the equivalent space-time code will 
not be the complete n x n Xd but instead will be missing 
some rows and will be of dimension (|D(5)| + 1) x n. For 
this, we now move to establish another property, necessary 
for the network’s approximate universality. We will name this 
property as ‘residual approximate universality’. 

Definition 1: Residual approximate universality is the prop¬ 
erty of an n x n approximately universal space-time code, 
which guarantees that after removing an arbitrary number, 
say k, of rows from each of the code matrices, it is then 
the case that the resulting ( n — k) x n truncated code is still 
approximately universal over any (n — k) x n r channel, for 
all n r . 

The proof that CDA-perfect codes are residually D-MG op¬ 
timal (and hence residually approximately universal) can be 
found in the proof of [11, Theorem 4] and is based on the 
fact that, given any truncated (n — k) x n codematrix X, the 
Hermitian nature of XX\ guarantees that the magnitude of 
each of its n— k ordered eigenvalues, is each lower bounded 
by the squared magnitude of the corresponding n — k smallest 
eigenvalues of the original pre-truncated matrix. 


A similar argument, gives that the above diagonal restricted 
perfect code is also residually approximately universal, over 
any MISO channel with i.i.d. fading. Having constructed a 
proper pair of residually approximately universal codes for 
the two stages, we can state that: 

Theorem 1: Given statistically symmetric Rayleigh fading, 
the optimal half-duplex constrained D-MG performance of the 
non-dynamic selection-decode-and-forward wireless network 
with 7i+l single-antenna users, is given by: 

rfND-SDAFopt(r) = (n - 1)(1 - 2r) + + (1 - r). (15) 

This performance can be achieved by utilizing a 1 x n 
horizontally-restricted perfect code during the first stage, hav¬ 
ing the intermediate relays decode if and only if they are not 
in outage with respect to the original transmitter, and finally 
having them re-transmit utilizing, together with the source, a 
distributed n x n diagonal-restricted perfect code. 

Proof: See Appendix ll .AI ■ 



Fig. 4. Optimal D-MG performance for the single-intermediate relay, ND- 
SDAF scheme utilizing perfect codes. Given a rate R (bpncu), cooperation 
applies for SNR > 2 2R . 


For results relating to different channel statistics and differ¬ 
ent channel topologies, we refer the reader to Appendix ll-BI 

III. Optimality in the Non-Dynamic 
Receive-and-Forward Strategy 

In this network, we initially consider the case where each 
node has a single receive-transmit antenna operating in half¬ 
duplex. The exact network introduced in Section II-C.2I devi¬ 
ates from the relay setup in [1], [2], as it requires one less 
intermediate relay since the destination D does not discard 
the direct transmission of S. Unlike in the case of the ND- 
SDAF strategy, the lack of information extraction at the 
relays, introduces the need that the related space-time code 
be completely distributable. Code distribution at the relays is 
discussed in Appendix II V- Al 

Interest in this ‘linear-processing’ (ND-RAF) relay network 
is generated by the scheme’s two main advantages. The first 
advantage is that the intermediate relays do not require channel 
information and the second advantage is that the intermediate 
relays do not perform time and energy consuming decoding. 
Potential power violation issues arising from utilizing the 
above properties are addressed in Sections GH and EH 










We proceed to further describe the relay network model and 
to analyze the equivalent second stage ‘two-product’ channel. 
The main result will be presented in Theorem |3 


A. Distributed space-time codes and the equivalent channel 

1) Relay scheme: During the first stage, the source’s single 
antenna sequentially transmits a vector 

k = 0z = 9-\zW ■■■ zW | (16) 

of n signals, where z is a codeword from a 1 x n coding 
scheme. Each intermediate relay It,. i = 2,3, ••• ,n then 
receives the n-length vector 

r. i = 0g i z + v i , (17) 

independently performs linear-processing (time-averaging 

based on space-time codes) on r ? ; and transmits 

= U A, (18) 

where each n x n matrix A,; is unitary. 

The signal at the receiver of the final destination D is then 
of the form 

n 

y = 'Y^hiX.i+w. (19) 

i= 1 

Clarifications regarding h\ can be found in Anncndix lll-AI 

2) Equivalent channel model: The authors in [1], [2] utilize 
a finite-length, random, distributed space-time code by having 
the linear-processing at the intermediate relays be performed 
by randomly chosen unitary linear-dispersion matrices [18], 
where each such matrix uniquely defines a codematrix row. 
The rationale behind this becomes clearer after rewriting < 1 1 91 
as: 

n 

y = y, hi(6gjZ + VjjAj + w 

i= 1 

n n 

= 0 hiQizAi + hi^Ai + w ( 20 ) 

2=1 2=1 

which easily transforms to the familiar point-to-point channel 
model 


where 


H 


y = 6HX + W 


( 21 ) 


gihi g 2 h 2 


Qnhn 


X = 


z_A\ 

zA 2 

zA n 


( 22 ) 


W = y hjVjAj + w. 

i =1 

Based on the results in Appendix EESI and for the sake of 
simplicity, we name the above channel as the ‘two-product 
channel’ by which we will exactly mean: 

Definition 2: The ‘two-product channel’ is the n x 1 MISO 
channel, modelled as in CD, with the channel matrix H as 
in CD representing fading coefficients which are products of 
two i.i.d. CA/"(0, SNR°) random variables. Furthermore, the 
effective additive noise consists of spatially and temporally 
white CA/"(0,1) random variables. 


B. Outage probability of the two-product channel 

Due to the absence of decoding at the intermediate relays, 
the rate-reliability limitations are mainly due to the relays-to- 
destination stage of the network transmission. As a result, the 
optimal D-MG tradeoff of the network is a function of the 
half-duplex effect and of the optimal D-MG tradeoff of the 
second stage two-product channel. 

The general method for establishing the D-MG limits of the 
two-product channel, relates to utilizing the existing eigen¬ 
value bounds of the approximately universal CDA codes, 
finding the codes’ ‘error contribution region’ (channel region 
in which perfect codes decode erroneously) and equating this 
region to the outage region of the channel. Establishing the 
eigenvalue statistics of the two-product channel then provides 
for the volume of the outage region and for the optimality 
limits. With respect to these limits, we note that: 

Remark 1: Since knowledge of g, at relay II, essentially 
reduces the two-product channel to the Rayleigh fading chan¬ 
nel, it is the case that the two-product channel’s optimal D-MG 
tradeoff cannot be any better than that of the Rayleigh fading 
channel. 

We proceed with the probability of outage for the two-product 
channel. 

Proposition 1: The optimal diversity-multiplexing gain 
tradeoff of the nxl two-product channel is given by 


d eq (r x ) = n{ 1 - r x ). _ 

Proof: From Lemma Bl (Appendix III-DL we see that for 

n 

A n =HH^=y WKfUf := SNR~ M 

i= 1 

being the only non-zero eigenvalue of the two-product channel, 
then the corresponding outage region, which at high SNR 
equals the error contribution region of the perfect codes, is 
given by: 

B = {p>l-r x }. (23) 

By definition of approximate universality, the probability 
of error outside B can be considered to be arbitrarily small, 
allowing us to limit our attention only to the channels with 
A n < (SNR 0 ) 2 . As a result, knowledge of the pdf of the 
channel outside this region is unnecessary, with the only 
condition that fx n (X n ) < oo, VA n >SNR°. We will see later 
in Annendix Ill-El that this condition is met. Immediately from 
Lemma |6] in Appendix EEH we see that for hi,gi, i = 
1,2,- ,n being CA/"(0,SNR°) random variables, then the 
probability density function of A„ = X]"=i ll^i|| 2 ||Si|| 2 ' s 
upper bounded as 

/a(A)<A" -1 (24) 


and for A„ = SNR~ M 

fM< SNR"' 1 ". (25) 


We here note that in [6], [9], the pdf of p for the nxl 
Rayleigh fading channel, is given by 

/^(/ x )=SNR~ Am e SNR ~ fI 

2 This is immediate by first considering that for the Rayleigh fading case 
r X)m ax = 1, and then by considering Remark Q 




which reduces to the two-product pdf expression of 
/m(/x)<SNR-^ for all p > 1 — r x > 0, that is for all p, € B. 
We then proceed as in [9], where B and J)Jp) completely 
defined the probability of codeword error as 

Pe< f f,Ap)d/i 

H&B 

since the double exponential nature of the probability of error, 
given a channel realization, acts as a binary indicator function 
in and out of B = {p > 1 — r x }. As a result 

P e < J SNR ~^ n dp 

MSB 

and using Varadhan’s Lemma [29] or the dominant term 
approach of Appendix II of [9], we get that 

P e < max{SNR _/i "}. (26) 

For a given multiplexing gain, the maximizing eigenvalue is 
then p = 1 — r x , which implies that 

P e < SNR~ n(1 - r * ) 

This lower bounds the optimal D-MG tradeoff in the nxl, 
equivalent two-product channel as 

d eq {r x ) > n(l — r x ). 

Remark [0 completes the proof. ■ 

We are now in position to give an explicit description of a 
D-MG optimal coding method for the ND-RAF network where 
the statistically symmetric Rayleigh fading channel is known 
only at the receiver of the final destination D, and where each 
node has a single transmit-receive antenna operating under the 
half-duplex constraint. The rate and the power constraint are 
known at all nodes. 

Theorem 2: (D-MG optimality in the ND-AAF and ND- 
RAF): The optimal half-duplex constrained D-MG perfor¬ 
mance of the ND-RAF network with n + 1 single-antenna 
users, is given by: 

dND-RAFopt(r’) = (n - 1)(1 - 2 r) + + (1 - r). (27) 

The performance is achieved by utilizing an approximately 
universal lxn horizontally-restricted perfect code during the 
first stage t = 1, 2, • • • ,n. If r > I then the relays do not 
forward the message and the source begins with the new 
message at f = n + 1. If r < i then each intermediate 
relay Ri, i = 2, ■ ■ ■ , n forwards at time t = i + n — 1 what 
it received at time t = i. Decoding at the final destination 
uses the received signals at either time slots t = 1,2, - • ■ ,n 
(r > ^), or at time slots t = 1 , n, n + 1 , ■ ■ • , 2n — 1 (r < ^). 

Proof: For r > the result is immediate by observ¬ 
ing that the equivalent channel is a SISO channel and the 
equivalent code is the approximately universal horizontally- 
restricted perfect code ®. For r < b the equivalent channel 
is the two-product channel, and the equivalent space-time 
code is the n x n diagonal-restricted perfect code d 1 41 that 
is D-MG optimal over all statistically symmetric channels. 
Consequently the proof is immediate through Proposition Q 


and from the fact that transmission takes place during 2n — 1 
time slots, shown in the proof of Theorem^to be the minimum 
allowed. ■ 



Fig. 5. D-MG optimal performance for the single intermediate relay, ND- 
RAF scheme vs the PEP based D-MG bounds for random and optimal coding. 


We now proceed to provide the general coding requirements 
for optimality in the ND-AAF and ND-RAF schemes, and to 
give bounds on the performance of networks that utilize some 
other existing coding methods. 


C. General coding optimality conditions in the ND-RAF relay 
network 


A closer look at the structure of standard perfect codes 
carrying n 2 information elements, reveals that approximate 
universality is a necessary but not a sufficient condition for net¬ 
work optimality, and that these standard-perfect codes do not 
provide for ND-RAF optimality as they fail to address the extra 
restrictions introduced by the half-duplex constraint. These 
extra constraints, exposed here in order to provide guidelines 
for future encoding methods, are a direct consequence of the 
fact that one cannot extract information at the intermediate 
relays, and can manipulate the received signals only with linear 
transformations. 

a) Condition 1: The distributed code must have at least 
one row that is the ordered subset of the vector transmitted 
by the source: This will allow for the destination’s decoder to 
utilize the signal of S. 

b) Condition 2: The n X T equivalent code should map 
exactly T information symbols from a discrete alphabet: 
Specific to the minimum delay case, we see in Lemma 0 
of Appendix IH-HI that annxT space-time code, carrying 
mT information elements from a discrete constellation and 
operating (in the point-to-point sense) at multiplexing gain 
r x , will allow for the network to operate at multiplexing gain 

On the other end of the spectrum. Proposition [8] tells 
us that the same code operating in the two-product channel 
with n-transmit antennas (m < 1), will achieve maximum 
multiplexing gain of r Xj max < rn. Consequently, optimal 
performance requires m = 1. 

Standard perfect codes, with m = n, fail this condition and 
provide for 


t^network-full (f) — {jl 1) ( 1 


+ n — 1 


+ (1 - r) 


n — 1 








guaranteeing for commencement of cooperation, given some 
rate R, only after 

SNR COO p.f u ii > 2 "- 1 R 

making it so that each new user entry renders the network less 
cooperative. 

c) Condition 3: Code does not require complex- 
conjucacy: The fact that complex-conjucacy is not a matrix 
operation, forces the source to send twice as many real 
information elements, thus doubling the duration of the first 
stage. 

For the single intermediate relay case, the Alamouti code, 
which can be shown to be optimal over the two-product 
channel, only provides for 

^network-alam() = (1 5 r)~ -\- (1 TJ 


guaranteeing for commencement of cooperation only after 


SNR, 


coop-alam _ 


> 2 


5 R 


in contrast to the optimal case, given by the diagonal- 
restriction perfect code of SNR coop . opt > 2 2R . 

For high values of n, orthogonal designs are expected to 
carry on the average rn ~ 4 information symbols per channel 
use and be of dimension n x n 2 . Analysis similar as that 
presented in the proof of Lemma[5]in Appendix III-DI (also see 
[23]), gives that the D-MG performance in the second stage 
is given by d{r) = n{ 1 — —) w n( 1 — 2r). The first stage 
duration needs to be of length 2 n 2 m ~ n 2 since conjugation 
cannot be described as a matrix transformation, providing for 

t^network-ortho (t*) ~ {jt 1) (1 T'f'j T (1 r), 


guaranteeing for commencement of cooperation only after 

SNR coop . ortho > 2 7R . 


D. Coding method for a network with a large number of users 
- the integral restriction perfect code 

We here present the ‘ integral-restriction perfect code ’ which 
manages to exhibit excellent performance especially as the 
number of users increases. Given a discrete information set 
A, the code is given by: 

XiT = {X = Efelo 1 /fc rfe > v /fc e A} 

where T from d60l provides the linear-dispersion matrices 
A,; = T n ~ l , i = 1, • ■ • , n. In addition to satisfying the 
necessary conditions of: 

1) raw data in one row 

2) one discrete information symbol per channel use 

3) no complex conjucacy 

it also manages to have: 

4) Constellation complexity remains constant as the num¬ 
ber of users increases: Due to the special nature of 
r (see [8], [10]) and after a small modification, the 
transmitted signals essentially belong in a normalized 


QAM-HEX signalling set, independent of the number 
of users. 

5) Fast encoding at the intermediate relays: Due to the 
sparse and discrete nature of the Af s, the scheme 
only requires one multiplication with a small Gaussian 
integer, per channel use. 

6) Network topology translates to reduction of the sphere¬ 
decoding complexity at the destination: The existence 
of only one receive antenna at the destination translates 
into a sphere-decoding complexity reduction at the final 
destination receiver, from 0{n 2 ) to 0{n). 

7) Allows for optimal rate, ease of construction and mini¬ 
mum delay, for any n. 

E. Communication with a base-station: reducing the half¬ 
duplex effect 

We will see that an increase in the number of antennas at the 
destination (base-station) does not only increase the diversity 
gain, given cooperation, but it also reduces the SNR required 
for cooperation to commence. This is done at the cost of extra 
decoding complexity at the final destination but bares no cost 
for the relays. 

In more detail, we recall a practical relay-network scenario 
given in [4] that talks of several relays with one trans¬ 
mit/receive antenna, cooperating in their task to communicate 
with a single base-station. It is logical to assume that the 
centrality of such a base-station will allow it to utilize multiple 
receive antennas. This can correspond to a wireless telephony 
setup where each mobile user utilizes the surrounding users to 
increase the reliability of the transmission to the base-station. 

It should be noted that the fact that the relays can only 
perform linear processing, prohibits having multiple antennas 
at the source since one cannot linearly process matrices in 
a meaningful way due to the additive nature of the received 
signal at the intermediate relays. Furthermore the intermediate 
relays can only have one receive-transmit antenna due to the 
lack of source-to-relay channel information. 

The following theorem explains how, given some network 
rate R, having multiple receive antennas at the base station 
can allow for a relay network with single-antenna intermediate 
relays to potentially reduce the SNR required for cooperation 
to apply. 

Proposition 2: Consider a base-station centered, ND-RAF 
relay network of n users cooperating through annxn space- 
time code whose D-MG performance over the equivalent 
(second-stage) channel is d eq (r). Let the code map on the 
average mn information elements. If the base station utilizes 
to receive antennas, 1 < rn < n, then cooperation is beneficial 
for multiplexing gains that are smaller than the multiplexing 
gain at the intersection of curves d eq (r(m + 1)) and d(r) = 
to( 1 — r). 

Proof: Since the information constellation satisfies 

I^I^SNR"* and since the entire network transmission has 
duration mn + n time slots, then under forced cooperation 
we have that 

-^network mn-\-n log 2 (SNR^) r% 

r ~log 2 (SNR) _ log 2 (SNR) _ to + 1 







Example 1: For n = 4 users cooperating to communicate 
with a base-station, doubling the number of antennas at the 
base station from m = 1 to m = 2 will allow for the necessary 
SNR, for cooperation to apply, to be reduced from SNR > 2 2 * 
to SNR > 2s*. 

The above bound hints towards utilizing n x n codes that 
map mn information elements and maintain sufficiently good 
eigenvalue bounds for increasing spectral efficiency. For this 
we turn again to the general family of CDA/perfect codes and 
consider the ‘m-layered n x n perfect code’ X m , given by: 

= {x = E^To 1 & (diag^ ■ G')) | 

which maps the mn information symbols from 
{/qi/E' - ’ i The codes have not been proven 

to be approximately universal. 

We conclude that one could accept an increase in decoding 
complexity and equipment, both only at the base station, in 
order to save power at the intermediate relays and to increase 
the SNR range in which cooperation is meaningful. 

We now move to the dynamic receive-and-forward cooper¬ 
ative diversity strategy. 

IV. Optimality in the Dynamic 
Receive-and-Forward Strategy 
A. Describing the scheme’s model 

For completeness, we will here reproduce the description 
and outage analysis of the dynamic amplify-and-forward net¬ 
work, first presented in [5] for infinite time duration, and we 
will then proceed to explicitly achieve this optimality in finite 
time duration. 

According to this strategy, the original source S transmits 
at each time instance, and each intermediate relay takes turns 
in transmitting an amplified version of a previously received 
signal. During a 2(n — l)-length frame, all intermediate relays 
have contributed, and the frame is repeated infinite times. The 
set of equations that describe each frame, as given in [5], is 

yt = hiXt+Wt, t - odd 

i = - + 1 

yt = hhi(giX i + + hix t +w t , t _ 2 yen 

or equivalently 

yt = hiXt+wt , f-odd 

Vt = bih igi x i + hix t + bihiV i}l + w t , 

where bi is the amplification factor at relay II,. For D-RAF, 
we set bi = 1, which, given a minimum and usually very small 
SNR, will always result in reduced average power consump¬ 
tion and thus no power violation occurs at the cooperating 
relays. 

Proposition 3: [5] Given infinite time duration, the optimal 
D-MG performance of the D-AAF and D-RAF schemes with 
n — 1 intermediate relays, is given by: 

do-RAF = (n - 1)(1 - 2 r) + + (1 - r). 

Proof: For completeness, the proof in [5] is briefly 
reproduced in Appendix IIII- Al ■ 


B. Optimal explicit construction for the dynamic receive-and- 
forward network 

We will here present an explicit construction, based on 
vectorized perfect codes, which achieves the optimal D-MG 
performance of the dynamic-receive-and-forward scheme and 
does so by using 2(n — 1) frames, each of duration 2(n — 1). 

Theorem 3: We consider the scheme where the source 
transmits continuously from a set xi, X 2 , • ■ • , £ 4 ( n -i) 2 
coming from a column-by-column vectorization of a 
2(n — 1) x 2(n — 1) D-MG optimal CDA space-time code, and 
where each intermediate relay It, , i = 2, • • • , n forwards at 
time t = 2 (i — 1) + 2(n — T)k, k = 0,1, • • • , 2 (n — 1) — 1 
what it received at time t = 1 + 2(n — l)/c. This n — 1 
intermediate relay D-RAF scheme achieves the optimal 
D-MG performance, as described in [5], of 

efo-RAF(V) = (n- 1)(1 - 2r) + + (1 - r). 

Proof: See Appendix IIII-BI ■ 

The above construction is optimal for any number n and 
for any statistical asymmetry. We here note that in [12], a 
CDA-based D-RAF scheme is proposed that, given statistical 
symmetry, achieves the optimal tradeoff with delay T = 
4 (n — 1). 

1) The emerging need for high-dimensional perfect codes: 
We have seen that in non-dynamic cooperative schemes, oper- 

R 

ating below some SNR threshold SNR coop = 2 rco °p will prompt 
the network to instruct the relays to not cooperate. If one 
insists on relaying even if SNR < SNR coop , then through CDA- 
based D-RAF, the same non-cooperative D-MG performance 
is expected (at higher decoding and signalling complexity). 
At low rates and high multiplexing gains, the corresponding 
SNR is bound to be small, rendering D-MG inaccurate and 
bringing us closer to outage capacity which relates to low SNR 
and which accentuates the role of maximization of mutual 
information. For a fixed rate, as the SNR increases, approxi¬ 
mate universality becomes more accurate and will guarantee 
for the fidelity of the results in TheoremsQtoQ] This high-low 
SNR duality emphasizes the need for approximately universal 
codes that are also information lossless, offer Gaussian-like 
signalling and tend to maximize mutual information, i.e. 
perfect space-time codes. This dual requirement, together with 
the fact that in D-RAF the dimensionality of the corresponding 
(28) codes grows fast with the number of users, jointly offer a 
substantial reason of existence for high-dimensional perfect 
codes [8], [10]. 

V. Comparing the presented cooperative diversity 

STRATEGIES 

A direct consequence of Theorems and0 is that: 
Theorem 4: The non-dynamic selection-decode-and- 
forward, the non-dynamic receive-and-forward, the non¬ 
dynamic amplify-and-forward, the dynamic amplify-and- 
forward and the dynamic receive-and-forward cooperative 
diversity strategies, provide for the same D-MG optimal 
performance. For all the above cooperative strategies, a 
network of n + 1 users, each having knowledge of the rate 




and power constraint, and each having a single transmit- 
receive antenna operating in half-duplex over Rayleigh 
fading, has a high-SNR probability of outage (optimal D-MG 
performance) P out := SNR ' iop, ' r / where: 

doptO) = (n - 1)(1 - 2r) + + (1 - r). (29) 

Optimality is achieved in minimum delay. 



Fig. 6. D-MG performance of the 3-node network for the above cooperative 
diversity strategies (ND-SDAF, ND-AAF, ND-RAF, D-AAF, D-RAF). 


VI. Simulations 

In Figure Q we present a comparison of the performance of 
all mentioned cooperative diversity schemes. The ND-AAF 



Fig. 7. Performance comparison for the 3-node (n = 2) D-MG optimal ND- 
SDAF, ND-RAF, D-RAF cooperative diversity schemes. The schemes are also 
compared with the 9-node ND-SDAF and the non-cooperative case. Network 
operates at 2 and 2.33 nbpcu. 


scheme, in which the amplification factor is based on the 
channel and set to equate the power consumption over the 
different frequencies, constantly performs worse than all other 
mentioned cooperation schemes of equal dimension. 

The plots assume that the users are asked to cooperate at 
any SNR. For the D-RAF, the amplification factor 62 was 
set to 62 = 1 , with the condition that the transmission over 
the second frequency (second row) did not violate the power 
constraint. This is always the case for the SNR range of 
comparison. The SNR recorded is the SNR representing the 
total power used, given that b 2 = 1. The same is valid for 
the ND-RAF case. The code used in the D-RAF case was the 
vectorized 2x2 perfect code operating over 4-QAM (M 2 = 4). 


In the n = 2 ND-SDAF, the code used was the 2x3 version 
of the diagonal restricted code, operating over 16-QAM, and 
in the n = 2 ND-RAF, the code used was the 2x4 version 
of the diagonal restricted code, again operating over 16-QAM. 
We note that this choice is not optimal but was used for reasons 
of network-rate uniformity. No information was allowed to be 
wasted at the receiver. In the n = 8 ND-RAF, the code used 
was the 8 x 8 integral-restriction perfect code. The decoding 
strategy of the first stage of the ND-SDAF network was such 
that the first stage horizontally-restricted code did not limit 
the error performance of the scheme. This was ensured by 
allowing the relay to decode only if 


1.921 2 > <5 


M 2 - 1 

SNR 


where 5 is chosen so that the effective rate is not the maximum 
allowed 8/3 bpncu, but instead 7/3 = 2.33 bpncu. This small 
drop in rate allowed for decoding to take place well out 
of outage, thus guaranteeing that the probability of error at 
the relay (given a horizontally-restricted perfect code and a 
given 5) is less than the probability of error at the receiver 
of the destination, given that cooperation took place. The 
choice of 5 and of the corresponding effective rates were 
arbitrary and somewhat heuristic. The following figure gives 
us an idea of the different choices and their utility. What we 



Fig. 8. Different decoding strategies at the intermediate relay, and their effect 
as to how fast cooperation is allowed to be fruitful. 


can conclude is that if the users have the choice to accept 
or reject cooperation then the preferred decoding strategy 
should be the one that allows for the maximum effective 
rate and intersects the performance of the cooperation scheme 
before the performance of the cooperation scheme intersects 
the performance of the non-cooperative case. On the other 
hand, if the users are forced to cooperate independent of 
SNR and rate, then the preferred decoding strategy presents a 
tradeoff between the effective rate and the SNR required for 
the probability of error at the intermediate relay to become 
smaller than the probability of error given cooperation. 

VII. Conclusion 

In this work we explicitly provided the hist ever D-MG 
optimal encoding schemes for several cooperative diversity 
strategies for wireless relay networks. These practical perfect- 
code based schemes are optimal over a broad range of channel 






















statistics and topologies. Bringing the results together we were 
able to conclude that the ND-SDAF, the ND-RAF (ND-AAF) 
and the D-RAF (D-AAF) have exactly the same D-MG optimal 
tradeoff. 

We proceed to recap the findings that relate to the different 
schemes. 

• For statistically symmetric ND-SDAF, ND-AAF and 
ND-RAF, optimality is provided by combining the 
horizontally-restricted perfect code over the first network 
stage, and the diagonally restricted perfect code over the 
second stage. This optimality holds for any i.i.d. fading, 
any number of users and with the least possible delay of 
T = 1, T = n or T = 2n — 1. 

• Optimality in D-AAF and D-RAF is achieved by 
letting the source transmit continuously from a set 
Xi,X 2 ,-'' , x 4 (n-i) 2 coming from a column-by- 
column vectorization of a 2 (n — 1 ) x 2 (n — 1 ) 
D-MG optimal CDA space-time code, and where each 
intermediate relay /?,, i = 2, ••• ,n forwards at time 
t = 2(i - 1 ) + 2 (n - 1 )k, k = 0, 1 , • • • , 2(n - 1 ) - 1 
what it received at time f=l + 2(n — l)fc. Optimality 
holds for any number of users and any statistical 
distribution of the fading. 

• For D-MG optimality, in a finite time duration ND- 
SDAF network, the first stage code needs to satisfy 
P(err | no outage) << P(outage), and the second 
stage code needs to maintain (the corresponding) D-MG 
optimality even when any number of its rows are deleted. 

- It was noted that D-MG optimal random Gaussian 
codes do not satisfy the first stage condition unless 
they are of infinite length. 

In comparing the ND-RAF, in a statistically symmetric 
network, with the other schemes we see that it : 


full-duplex D-MG performance will be optimal over any 
channel with i.i.d fading, and will be given as, 

P en = P Xh (l,l)(y / P(\D(r 0 )\=i)P Xd (i,l)\ 
o ' 

= Pa* t(l, 1) (]T P(|PM| = t)Pout(i, 1)) 

even for the minimum delay case, where P out (fc, l) represents 
the probability of outage in the k x l channel. For the Rayleigh 
fading MISO channel, we have that 

Pout(», 1) = (Pou.(l, 1))' = SNR- i(1 " r) 

resulting in the full-duplex probability of network transmission 
error, for the statistically symmetric Rayleigh fading case, to 
be given by 

Pen = P 0 ut(l,l)(^(Pout(l,l)) n “ 1 “ i (Pout(l,l)r) 
M=0 ' 

= (Pout(l,l)) n = SNR-^ 1 -^. 

In regards to the half-duplex constrained D-MG performance 
in the statistically symmetric case, we observe that the overall 
network transmission duration was 2n — 1 time slots. The task 
is now to prove that the duration cannot be any less. For this 
we first turn to Proposition [8] which tells us that to achieve a 
maximum (second stage - MISO) multiplexing gain of r max = 
1, we need to map n discrete symbols in A’,/. For this we need 
at least n time slots in the first stage (SISO) since any smaller 
duration would violate the r max = 1 limit placed by the first 
stage SISO channel. The proof is then complete by observing 
that we need at least one time slot for every row of Xd that 
does not correspond to direct transmission. 


• requires no channel knowledge at the intermediate relays 
(compared to full knowledge for the ND-SDAF) 

• has minimum delay of 1 or 2n — 1 (compared to T = 
4 (n — 1) of the least-delay D-AAF scheme [12]) 

• has signalling complexity of n (compared to complexity 
of 2(n — 1) for the D-AAF scheme) 

• has sphere decoding complexity of n (compared to 
4 (n — 1) of the best proposed D-AAF scheme). 

As a final thought, it is pointed out that in terms of D-MG, 
what was believed to require complete channel knowledge 
at the receivers of the intermediate relays, require encoding 
over infinite time duration and require decoding of infinite 
complexity, was exactly achieved by newly constructed perfect 
code variants, given absolutely no channel knowledge, given 
minimum time duration and small decoding complexity. 


B. D-MG tradeoff for different fading distributions and differ¬ 
ent channel topologies 

a) Optimal D-MG performance for some other fading 
distributions: Due to statistical symmetry, it will be the case 
that, for any channel probability density function (pdf), the 
network’s D-MG optimality requires the same outage-based 
decoding strategy 0 , as well as the joint element property 
and the residual approximate universality of the utilized space- 
time codes. For an exact network performance expression, 
we also need to know the codes’ D-MG performance over 
the corresponding channels. Deviating from Rayleigh fading, 
analysis in [7, Prop 5.2] provides the optimal tradeoff curve for 
the class of n x 1 MISO channels with i.i.d fading coefficients 
{cj” =1 , to be 


Appendix I 

Results and proofs relating to the ND-SDAF 

SCHEME 

A. Proof of Theorem 0 

Given the approximately universal Xh and the residually 
approximately universal Xd, it is the case that the network’s 


d*(r x ) = an(l — r x ), 


0 < r x < 1 


where a := lim los ( p ,^ c ,^ —We conclude that the overall 

t—o lo s( 4 ) 

network half-duplex D-MG optimal performance in such a 
channel is 


(WsDAFoptM = a(n - 1)(1 - 2 r) + + a(l - r). 



1) Encoding schemes for networks with nodes having mul¬ 
tiple receive and transmit antennas: We consider the special 
case where each terminal has to receive-transmit antennas, 
and where the fading across the different antennas has an 
arbitrary statistical distribution , with this distribution being 
the same over any path. This topology brings to the fore the 
’horizontally-stacked perfect code’ and the following bound. 

Proposition 4: The full-duplex D-MG performance of the 
ND-SDAF with n + 1 nodes, each having m transmit-receive 
antennas, is given by the integral-point-wise linear plot gov¬ 
erned by: 

d*(r x ) = n(m - r x ) 2 , r Xjmax = to 

and the half-duplex performance is lower bounded by the 
integral-point-wise linear plot governed by: 

rfND-SDAFm(r) > n{m - 2 r) 2 , 0 < r < ^coop 

= (to - r) 2 , r C oop < r < TO. 

where r coop is the intersection of integral-point-wise linear 
plots n(m — 2r) 2 and (to — r) 2 . The bound is met by 
utilizing an to x m(n — 1) horizontally-stacked perfect code 
(horizontal stacking of n — 1 independent matrices from an 
to x to perfect code) during the first stage, and a distributed 
m{n — 1) x m{n — 1) perfect code over the second stage. 

Proof: See Appendix ll-CI ■ 

The proof of approximate universality of the horizontally- 
stacked CDA code is found in [11, Section V.B], and is based 
on the fact that the Hermitian nature of XX ' guarantees that 
the magnitude of each of the to ordered eigenvalues of any 
XX\ is lower bounded by the magnitude of the corresponding 
eigenvalue of any of the X^ k \ k = 1 , 2, - - - ,n- 1 . 

2) Approximately universal codes and the most general 
relay network setup: We here note that, given any intermediate 
relay decoding strategy, it is the case that the generalization of 
the above multi-antenna, multi-node, CDA-based scheme, will 
allow for an overall network full-duplex D-MG performance, 
that is optimal over all other point-to-point encoding schemes, 
even for the most general case of having 

• an arbitrary number of receive antennas varying over each 
terminal 

• an arbitrary number ns and n ti of transmit antennas 
varying over the different nodes S and { R, }"_ 2 

• arbitrary and possibly different probability distribution 
functions for the fading across the different antenna pairs 

• an arbitrary number of users 

• an arbitrary decoding strategy 

• an arbitrary network transmission duration. 

Specifically, for some T depending on the network setup, 
this generalized optimal network encoding scheme asks for 
an ns x T CDA code to be used for transmission during the 
first stage, and which will be part of a {ns + Sr=2 n i») x 
{ns + 2 n ti ) CDA code that will be distributed across all 

relays, for transmission during the second stage. 


discrete constellation A m of cardinality A rn = SNR” 1 , 
into 7i—l codematrices {A'( fc )}^), each from an approx¬ 
imately universal CDA space-time code X m of cardinality 
X m = SNR rm . Having done so, S sequentially retransmits 
0[X« A'( 2 ) • • • X( n ~ 1 ' ) ], essentially sending a to x m{n — 
1) codematrix X = [X^ A'( 2 ) ••• X^ 11-1 )] from the 
horizontally-stacked CDA code whose approximate universal¬ 
ity was proved in [11, Section V.B]). The resulting guarantee 
that P{ err | no outage ) = SNR~°°, makes it such that, by 
t = (n — l)m, all intermediate relays in D{S) have cor¬ 
rectly decoded X and gained full knowledge of {/i}™ /” 
Consequently, each relay in D{S) now re-maps the correct 
information into {m{n — l)) 2 elements from a discrete con¬ 
stellation of cardinality = SNR™!”- 1 ), 

and eventually re-maps these new information elements into 
the codematrix X m r n _ i) of an m{n — 1) x m{n — 1) CDA 
code. 

In the second stage, the k th transmit antenna (k = 
1, • • • ,to) of intermediate relay Rj £ D{S), transmits the 
{{j — 1 )m + k) th row of codematrix X m ( ra _!). As a result, the 
participating relays essentially construct a (|£>(S , )| + 1)toxto7i 
punctured version of the mnxmn CDA code. Due to the resid¬ 
ual approximate universality of the CDA codes, we have that 
this (.D^Ito x mn punctured code is approximately universal 
over any \D{S)\m x n r channel, for all n r . Furthermore we 
know that the optimal decoding strategy will ask for the relays 
to decode if and only if the corresponding to x to channel is not 
in outage. This, combined with the fact that the codes, utilized 
on both stages, have the required approximate universality and 
residual D-MG optimality property, allows for the proposed 
encoding scheme to provide for the optimal full-duplex D-MG 
performance over the network. 

Narrowing to the specific cases of the Rayleigh fading 
channel, we first observe that for i = 0,1, ■ ■ ■ , n — 1, it is 
the case that the expression 


P(|D(S')| = i)P out {im, to) 

= (P 0 ut{m,m)) n 1 l P oat {im, to) 

• (m—r)(m—r)(n— 1 — i)— {inn— r)(m—r) 


is maximized for i = 0. Consequently, using the dominant 
summand approach as in [6], [9], we conclude that the optimal 
full-duplex D-MG performance over the network, is given by 

( n— 1 

^P(|D(5)| = i)P oat {im, to) 

2=0 

= P 0 ut{m,m)P{\D{S)\ = 0) 

= P 0 ut(TO,?n)(P 0 ut(TO,TO))" 1 
= (Pout(TO, TO))” 

= SNR~ n(m - r)2 . 


C. Proof of D-MG bound on the multi-antenna ND-SDAF 
scheme (Proposition 0) 

The encoding scheme first asks for the source S to 
map m 2 {n — 1) information elements {/»}?L^ n_1 \ from a 


The fact that 2 m{n— 1) time slots are required for the network 
to complete the transmission of the m 2 {n — 1) information 
symbols, from a constellation A rn of cardinality A rn = 
SNR”*, concludes the proof. □ 



Appendix II 

Proofs relating to the ND-RAF scheme 

A. The equivalent second stage ‘two-product’ channel 

We first note that we can equate having the first linear 
dispersion matrix Ai be equal to the identity matrix I n , to 
correspond to the direct uncoded transmission between the 
source and final destination. Decoding with consideration of 
the source’s signal, results in savings of one intermediate relay. 
Having set A\ = I n translates to considering the source as an 
intermediate relay, with g i = 1 known to the receiver and the 
transmitter. For uniformity we will assume that the transmitter 
does not know gq. As we proceed, it will become apparent 
that this does not affect our analysis. Consequently, we will 
henceforth consider gi,hi, i = 1, • • • ,n to be i.i.d, zero- 
mean, complex Gaussian random variables. As in [1], we now 
consider the channel equations E3 and 03, and observe that 
matrix X can be considered as an n x n codematrix and matrix 
H as the 1 x n channel fading coefficient matrix. The authors 
in [1], [2] show that due to the required unitary nature of the 
linear-dispersion matrices and due to the fact that the hi ’s are 
known to the receiver, it is the case that the effective additive 
noise elements W\hi are spatially and temporally white, zero 
mean, Gaussian random variables, with variance 

n 

Var{W\hi) = (l + J2 N 2 ) J »- (30) 

2=1 

In Appendix III-BI we see that for unitary Ai s, in the 
SNR scale of interest, the W\hi can be considered, without 
loss of generality, to be CA/"(0,1) random variables, con¬ 
ditioned on considering the gi,hi fading coefficients to be 
i.i.d. CA/"(0, SNR 0 ) random variables. This is exactly the two- 
product channel. □ 


B. Proof of the i.i.d nature of fading 

As stated in [1], [2], W\hi are spatially and temporally 
white CA/"(0, (1 + Y^i=x |^i| 2 )) random variables. The prob¬ 
ability density function of h = ]P’ l =1 |^z| 2 is 

h n 


p{h) = 


i-i e - h 


(n- 1)! ' 


For e > 0, 


p(h = SNR £ ) = 


SNR 


(n-rj^-SNR- 


(n — 1)! 

and due to the double exponential term we have that 
p(h = SNR £ )=SNR~°° = 0. On the other hand 


p{h = SNR" e ) = 
and since e 


SNR 


-(ra-l)e e -SNR 


(«-!)! 

-SNR we have that 

p{h = SNR~ £ )=SNR~ £(n_1) 


which implies that Har(VF|/ii)=SNR°. 

Now we observe that 

E[\\H9X\\ 2 F ] _ CKTD _ n r n t Var(gihi)¥.[\\9X\\p] 


n\\w\hi\\ 2 F ] 


= njSNR = 


n r TV ar{W\hi) 


and thus 

Finally we note that in the SNR scale of interest, we may 
interchange the SNR 0 and 1 values and have Var(gi) = 
Har(/ii)=SNR° and Var(W\hi) = 1, since the substitution 
maintains the power-SNR requirements. 


C. Existing performance bounds for the ND-RAF (linear- 
processing scheme) (Proposition |6j 

a) The diversity result from [1 ]: Having established the 
whiteness of the additive noise, the authors in [1], [2] proceed 
to prove the following: 

Proposition 5: (Variant of the main result in [1]) Over the 
effective channel described in (ED, with the channel fading 
coefficient matrix as in E3, utilizing full-rank matrices as in 
<E3 and with white additive noise as in <E3, it is then the case 
that in the high SNR regime, the pairwise error probability 
based diversity is n. 

In expanding the above PEP bound to consider multiplexing 
gains other than r = 0, we get: 

Proposition 6: Given random coding, the existing PEP 
bound in Proposition 0 translates to a D-MG bound (of the 
entire network) of 

((nd-raf(t) > (n - 1)(1 - kr) + + (1 - r), k » 1. 


Given coding based on approximately-universal codes, the 
same PEP bound translates to a D-MG bound for the entire 
network of 


g(nd-raf(t) > (n — 1) ( 1 — 


An — 1 
n — 1 


+ (l-r). 


Proof: It is not difficult to see that in the high-SNR 
regime, the PEP diversity of the second stage of the network 
scheme, as given in [1], [2], is a function only of SNR and 
of the minimum eigenvalue of any difference of any two 
codematrices (at a given rate), and is upper bounded as: 


PEP<SNR~ n (l' n )- n (31) 


where l' n is the smallest eigenvalue of the power normalized 
version of the code where the codematrices X' are of the form 
X' = ^ such that E[||X , |||,]=1. 

For the channel model given as Y = 6HX + W, 
with fE[||J 5 G||^ r ]=£' 2 and 6 such that E[||0X|||.]=SNR, 
it was shown in [9] that a non-vanishing determinant 
det(AXAX^)>SNR°, guarantees a minimum partial eigen¬ 
value product of 

n 

n u > (e 2 )-^- 1 ) 

i=n—j 

and minimum code-eigenvalue of l n > ( E 2)-("-D. 

Interestingly, in the normalized case of E(||A |||.) = 1 with 
f = jf;, the minimum partial eigenvalue product becomes 

n 

n ^>(^ 2 r n vj 

i=n—j 









which accentuates nicely the fact that code performance de¬ 
pends solely on its smallest eigenvalue. For when the n x n 
CDA code X maps n 2 information symbols from a discrete 
constellation A, we have that 

E 2 = \A\ = \X\& = (: 2 Rn )£ = ( 2 nrl °s SNR )^ = SNR-. 

As a result ^>(SNR") -n = SNR -r which brings us to 

PAP<SNR - "SNR -r 

The overall probability of error is then 

P e < \X\ (PEP) = SNR nr SNR - "SNR -r 

which allows for lower-bounding the D-MG tradeoff of the 
two-product channel as 

rfequiv(f) > n — 2rn (32) 

implying a maximum diversity d e quiv(0) = n and maximum 
multiplexing gain of r eq uiv > 5 . Consequently, a temporally 
disjoint first and second stage result in cooperation applying 
at the solution of n( 1 — 4r) = 1 — r. This concludes the proof. 

■ 

Given approximately universal codes, the PEP based D-MG 
bound can only guarantee that cooperation will apply for r < 
which, given some rate R, translates to having coopera- 

An — 1 p 

tion only after SNR > 2 "- 1 . Given finite duration random 

codes, cooperation might require SNR > 2 kR , k » 1 . 

D. Lemma on the error contribution region of perfect codes 

Lemma 5: For A„ = HH^ = ^=1 |N| 2 ||ft || 2 := SNR 
being the only non-zero eigenvalue of the two-product channel, 
then the corresponding outage region which at high SNR 
equals the error contribution region of the perfect codes, is 
given by: 

B = {p>l-~r}. (33) 

Proof: Directly from the mismatched eigenvalue theo¬ 
rem presented in [9], we know that the minimum Euclidean 
distance between any two codematrices after the action of the 
channel, is lower bounded as 

4 = \\0hax\\% > e 2 \j n (34) 


simply because we have a double exponentially decreasing 
term multiplied with a polynomially increasing term. As a 
result, we have the error contribution region, being defined as 
the set B of eigenvalues p such that P{err\p) fz 0. This 
region requires that c < 0. Consequently, gives that 
c = 1 — r — p < 0 which concludes the proof. ■ 

E. Lemma on the effective pdf in the two-product channel 

Lemma 6: For hi, gi, i = 1,2, • ■ • , n being CA/"(0, SNR 0 ) 
random variables, then the probability density function of 
A„ = Yh= 1 141| 2 1|&II 2 is upper bounded as 

fx( A)<A n " 1 (36) 

and for X n = SNR -/i 

fM< SNR-' 1 ". (37) 

Proof: We first recall that for x, y being two inde¬ 
pendent random variables and for z = xy, we have that 

fz,x{z, x) = i^y f x ,y(x, f) and thus 

OO 

fz{z)= [ ^-f x ,y{x,-)dx 

J M x 

— OO 

Now consider hi,gi being two independent identically dis¬ 
tributed CA/”(0, k), A:=SNR°, random variables and thus 
x = |4|| 2 , V = IISi|| 2 > two i-i-d. exponential random variables 

with f x (x) = e~ kx , f v (y) = e~ ky , f x , y (x,y) = e - kx e~ ky 
and f x , y (x, |) = e~ kx e~ k x, x,y,z € R + . As a result, 

OO OO 

f z {z) = [ -e~ kw e~ k ^dw = [ -e- k{w+ ^dw. 

J w J w 

0 0 

From equation 3.471.9 of [26] we see that 

f z (z) = 2 A' 0 ( 2 fc v / i) i.e. 

f\\hi\\ 2 \\gi\\ 2 ( z ) = 2 A' 0 ( 2 fc v / i) 

where Kq(-) corresponds to the modified Bessel function of 
the second kind. 

From [27, Section 6 . 6 ], we observe that 


where l n corresponds to the smallest code eigenvalue of 
AXAJ4 For the case of the perfect code X, which maps 
n 2 information symbols from a set A, we know that 0 2 = 
SNR 1- ", („>SNR“" (n_1) . As a result, 

d 2 E {p) > SNR 1 - "SNR~ n SNR _A1 := SNR 0 . (35) 

Observing that the pairwise error probability, given a channel, 
is upper bounded as 

PEP(p) < 

allows for the entire probability of error, given a channel p, 
to be upper bounded as: 

P(err\p)<\X\ ■ PEP(p) = SNR" r e" SNRC . 


It is clear that 

lim SNR" r e" SNRC 

SNR->oo 


Ko(x) 

lim ——— 

a —>0 —Log{x) 


1 . 


Observing that for x>0.5, K 0 (x)<^y=^e~ x < 1 is always 
finite (and decreasing in a double-exponential rate), allows 
one to interchange SNR 0 and 1. Using the absolute unity, we 
also note that as e —> 0 then 


K 0 (l - e) 
~log{ 1 - e) 


-1 


We then choose a 

Jfo(1 ~ e) = fc=SNR° 


guaranteed to exist e such that 
-iog(i-e) — 1 ji'i' ■ It is easy to see that 1 — e is closer to 
1 than to 0. This combined with the monotonically decreasing 
nature of I\q(x), gives that 


1 —e 


Ko{x)dx> / Ko(x)dx 


= SNR -00 =0, Vc> 0, 


0 


1 —e 







and as a result of the dominant summand effect 

1 1—e 


Ko(x)dx= j K o (x)dx + 0. 


Consequently, the term J^_ t Ko(x)dx can be substi¬ 
tuted by the smaller term _ — log(ir)(i:E. Combined 
with the fact that /\ 0 (1)=A' 0 (SNR 0 )=SNR 0 and that 


Kq(xi 

- i°gOi; 

we get 


< 


Kq(x 2 ) 
- log(:Z2) 


\/xi < x 2 , it is the case that for all z< 1 , 


z z z 

J I<o(\fz)dz< — k J log {y/z)dz= — J log (\fz)dz 


Recall that if A £ B then A<SNR°. The fact that all the 
summands of the channel eigenvalue A = |<7i| 2 |^i| 2 are 

positive, suggests that if A £ B then |g,;| 2 |/ii| 2 <SNR°. This 
justifies limiting our attention to the f\ g .\ 2 \ h .\ 2 (x), x < SNR 0 . 
As a result, for hi,gi being CA/"(0,SNR°) random variables, 
then the cumulative distribution function of ||/ii|| 2 ||<?i|| 2 , up to 
SNR 0 , is upper bounded as 

^(IIM 2 |bl| 2 <*)<z(l-l°g(*)) (38) 

since 

Z Z 

F z (z) = - J log(Vt)dt= - i J log(f)df 
0 0 

= i(f-flog(f))|g = z — z log(i:) + 0 log(O) 
and since 0 log(O) = 0, we get that 

p (IN| 2 ||5il| 2 < z)=z{ 1 - log(z)). (39) 

To complete the proof of the main Lemma [ 6 ] we use the 
following proposition: 

Proposition 7: Consider n independent identically 
distributed random variables ,a n with a,; > 0, Vi. 

It is then the case that for A = ^" =1 a * then 
P(\<t)<[P( ai ) <t] n . 

Proof: Sketch of proof of Proposition |7] In the n = 2 
case, the ( 01 , 02 ) pairs such that ai + 02 < t are the points 
inside the triangle formed by the positive ai-axis, the positive 
02 -axis and the line 02 = t — a\ = f(a\). As a result, for a 
certain a' 2 , it is the case that a’ 2 + a\ < t , Voi < f~ 1 (a 2 ). 
The corresponding triangle is completely enclosed in the t x t 
square whose lower-left corner is at the (0, 0) point. As a result 
P(oi + a 2 <t) < P(oi < t)P(a 2 < t) = (P(oi < i)) 2 . For 
the 3-dimensional case, the triangle is replaced by a pyramid 
which is enclosed in a 3-dimensional cube with edge of length 
t, and in the general n-dimensional case the n-dimensional 
pyramid is enclosed in an n-dimensional hypercube again with 
edge of length t. This concludes the proof of Proposition Q ■ 
Given the above proposition we continue from d39l and 
can now deduce that, for the cumulative distribution func¬ 
tion P(||/iiII 2 1 ||| 2 < z) < z( 1 — log(^)) presented in equa¬ 
tion (1381 and for A = IIM| 2 ||Vi|| 2 ’ we h ave that 

P(A < z) < [P(|N| 2 |N| 2 < *)]" < [z{l - log{z))] n 


Consequently 

Pa(A) < [A(l — log(X))] n 

Keeping in mind that the range of interest is 
A < SNR 0 and that P A (SNR°)=SNR°, then for 
A = SNR' 5 , SNR —> 00 , 1 > <5 > 0, it is the case that 
(SNR' 5 (1 - log(SNR‘ 5 ))) m = (SNR' 5 )™, m > 0. As a result, 
the cumulative distribution function (cdf) is expressed as a 
polynomial with only one term as 

Pa(A)<A™ 

which allows us, after differentiation, to apply the existing 
upper bound inequality of the cdf to the pdf, and get that 

/a (A) < A "" 1 (40) 

For A = SNR -/ \ p = then f^p) = ±4^ with 

t = -i^SNR) A' As a result ’ 

fM < SNR-^log(SNR) (41) 

As in [ 6 ], the fact that 

ita 1 ° g « 1 °f<™ R >» -»0 

SNR^oo log(SNR) 

indicates that the log(SNR) term does not contribute to the 
SNR exponent of the probability of error, and as a result 

fM < SNR"'*’* (42) 

which concludes the proof of Lemma [ 6 ] ■ 

F. Proof of the antenna limitations of the relay network 

The fact that the relays can only perform linear processing, 
prohibits multiple antennas at the source since one cannot 
linearly process matrices in a meaningful way due to the 
additive nature of the received signal at the intermediate 
relays. Furthermore the intermediate relays can only have one 
receive-transmit antenna due to lack of source-to-relay channel 
information. 

G. Relation of number of mapped symbols to r max 

Proposition 8: Consider a full-diversity n x T, T > n, 
space-time code X that maps, through linear combining, on the 
average mT information elements from a discrete constellation 
A. Given that the code operates in the Rayleigh fading with 
n-transmit and n r receive antennas, with n r > m, or operates 
in the two-product channel with n-transmit antennas (to < 1 ), 
then such a code can achieve a maximum multiplexing gain 
of r max < to. 

Proof: We have that |Tj = 2 RT = 2 rT j 0 ^SNR = 
SNR rT = |Al| m which implies that |Al|=SNR™ and since 
the constellation is discrete E[||a £ _4|| 2 ] = |Al|. The fact that 
each element X i 3 of a code matrix is a linear combination, 
with coefficients independent of SNR, of elements of A, 
gives us that IE[|| 1 | 2 ] = \A\ = SNR™. For 9 such 
that E[\\9HX\\ 2 f \ =’E[||0W|| 2 .]=SNR we have that 9 2 = 
SNR 1- ™. 








Without loss of generality we can assume that there exist 
two codeword matrices Xi , X 2 £ X, with X, mapping the 
information nm-tuple {a,;, 0,0, • • • , 0}, where ctj=SNR°. As 
a result, the determinant and trace of the difference matrix 
AX, is a polynomial of degree less than n over a = a± — 
a 2 =SNR°, with coefficients independent of SNR, i.e. 

det(AXAX t ) = Tr(AXAX t ) = SNR 0 

and thus with all its eigenvalues l, = SNR 0 . The corre¬ 
sponding pairwise error probability PEP(X —> X'), in the 
Rayleigh fading channel, then serves as a lower bound to the 
codeword error probability P e , i.e.. 


P e > PEP(X -> X ) = 


1 


n? =1 [ i + fy 


= SNR 


—n r n(l— - 


dir) < n r n (1 -) 

V 771/ 


(43) 


Appendix III 

Proofs relating to the D-RAF scheme 

From [5], we see that due to the independence of the 
information symbols, the power constraints are satisfied when 


bi < 


SNR 


IffipSNR + 1' 


We denote by <fii and A the exponential orders of h., and fj t 
respectively, i.e. 


hi= SNR"^, gi = SNR" 4 * 

and as in [5], we observe that 

,, , / SNR" 00 , &<0 

PM&i) | SNR - *, > 0 ' 


(45) 


(46) 


which proves the claim for the Rayleigh fading case. Remark 
[2 concludes the proof. ■ 


The same holds for Asa result, in accordance with [ 6 ], [9], 
[5] and without loss of generality we will limit our attention 
t° £ R + . Given that the probability that £ R _ 

can be considered to be arbitrarily small, we may restrict the 
amplification factors to be 6, = 1. 


H. Relation between rate reduction due to half-duplex and 
number of mapped symbols 

Lemma 7: Consider a distributed n x n space-time code, 
with one of its linear-dispersion matrices equal to the identity 
matrix, and with D-MG performance over the two-product 
channel given by d eq (r). Furthermore assume that the space- 
time code carries mn information elements from a discrete 
constellation. It is then the case that the D-MG tradeoff of the 
entire half-duplex constrained linear-processing relay network 
utilizing this code, is given by 


dnetwoi-k(r) = d eq (r(m + 1)). (44) 

Proof: We first observe that linear-processing does not 
allow for the intermediate relays to extract information from 
received signals since the distributed code X can only hold 
linear combinations of the received symbols. This forces the 
source to transmit the mT information symbols one at a time, 
and the intermediate relays to transmit over T time slots, 
without information extraction and re-encoding. This means 
that the lack of decoding at the intermediate relays does not 
translate to arbitrarily high rates. It is in fact the case that the 
rates are limited both by the capacity-outage corresponding to 
the dimensions of the second stage channel, as well as from 
the fact that the transmission during the second stage will have 
to take place over = T time slots. 

For r corresponding to the multiplexing gain of the second 
stage, it is the case that \X\ = SNR rT and as a result the 
information constellation satisfies 


A. Proof of Proposition 0 

Beginning with the single intermediate relay case, a frame 
(El}, is given by 

y=[yiV2] = s + n 

= [h\X\ h 2 g 2 xi + hix 2 \ + [wi h 2 V2,i +w 2 ]. 


For compactness we will denote SNR by p. For 


G 2 


hi 0 
h- 2 g 2 hi 


(47) 


denoting the equivalent channel defined by the above equa¬ 
tions, we have E s = pG 2 G! 2 denoting the covariance matrix 
of the signal and fading part of the observed vector at 
the destination, with the expectation taken over the signal¬ 
ing set and the fading coefficients held fixed. Furthermore, 


— 


denotes the covariance matrix of the 


1 0 
0 1 + \h 2 \ 

additive noise part of the observed vector at the destination. 
Using the fact that the largest eigenvalue Z max of E^a;*] 
satisfies l max =p and directly from El}-®, the mutual in¬ 
formation during a single 2-length frame is given in [5] to 
be 


I(x;y\hi,h 2 ,g 2 ) = log 2 det (/ 2 + pGG t E„ 1 ) (48) 

= max{ 2 (l-(/>i) + ,l-^ 2 -£ 2 }log 2 p. 


|yf|<SNR-. 

Consequently, in the T = n case, for d(r) corresponding 
to the transmission over the n time slots of the second-stage, 
it is then the case that d(r(m + 1 )) corresponds to the entire 
mn + n time-slots (first and second stage duration). ■ 


We now consider the outage event, that is the event where 
the corresponding maximum mutual information allowed by a 
set of fading coefficients is less than the rate of transmission 
R = r log 2 p. That is 

O = {(hi,h 2 ,g 2 ) £ C 3 | I{x;y\hi,h 2 ,g 2 ) < RT f } 









where Tj is the duration of the frame. Furthermore, due to and modify to reflect noise whitening by 
d 46 > . we are only interested in 


0 + = 


E = E 


= j(0 1,02,6) e K +3 | max{2(l-0i) + , 1-02-&}} < 2rj. 
The probability of outage is given by 

Pout(r) = p“ rf ” (r) = f PMufai&dfadfodb 


Wl 

h 2 V 2 ,l + W 2 
1 0 
0 \h 2 \ 2 + l 

the equivalent whitened model 


[ w t h : 


2 u 2,l 


Wo 


( 56 ) 


i-i 


y i 2/3 

2/2 2/4 


sup p 

(<h,<h,ta)€0+ 


(0i,02,€2)eo 

-( 01 + 02 + 52 ) 


( 50 ) 


The last equation comes from the dominant term approach in 
[ 28 ], extensively used in [6][ 5 ][ 9 ], or directly from Varadhan’s 
lemma [ 29 ]. Consequently 

d 0 ut(r) = inf {01+02+£2} ( 51 ) 

max{2(l — 0i)+,l — 02— 52}<2r 

which readily leads to the final expression 

cfc)-RAF(r) = (1 - r) + (1 - 2 r)+. 


= E" 

+ E ' 1 


1 

hi 0 


Xi x 3 


1 

<N 

(N 

C53 


X 2 X 4 


Wl 


w 3 


h 2 V 2 ,l + w 2 h 2 V 2>3 + W4 


( 57 ) 


For 

E" 

and 


1 

Wl 

w 3 


Ml M 3 


^2^2,1 + 

h 2 v 2t 3 + W4 


M 2 M4 


Now moving to the multiple intermediate relay case, we 
again assume Gaussian random coding across the frame, and we a i so ] e t 
introduce E Si , i = 2 , 3 , • • • ,n 

E Si = pGiG\ = p 


hi 

0 1 


r *1 

0 1 

_ 52*2 

hi 

— 

Q2 h 2 

hi 

ill 


L 1++2I 2 

I++2I 2 J 


l ft i I 


h\h l+ ig i+ i \hi \ 2 + \h i+ 1 \ 2 \g i+1 \ 2 


Zl 

z 3 

= E _1 

yi 

53 

. ^2 

Z4 _ 


. 52 

54 _ 


and E n( , i = 2 ,3 

En . 


*1 

0 

xi x 3 

Ag 2 h 2 

/11A 

x 4 




: 


0 


0 a 2 w + \h i+ 1 \ 2 a 2 0 l + \h l+1 \ 2 

As a result, the block diagonal covariance matrices are now where A = 


+ 


Ml M 3 
U 2 U4 


E* = 


Ss 2 02x2 
02 x 2 


j — 


02 x 2 

02 x 2 


( 52 ) 


1++2I 2 
Heir = 


. We now define 


/ii 0 

Ag 2 h 2 h\A 


Proceeding as in the single intermediate relay case, concludes 
the proof for the probability of outage corresponding to the 
equivalent channel 


G = 


G 2 02x2 

02 x 2 G 3 


42x2 


0 2 


X 2 


02 x 2 

02 x 2 

G n 


and proceed to find the outage probability by calculating 

det(J + = 

l + p\hi \ 2 phi{Ag 2 h 2 )* 

p{Ah 2 g 2 )h\ 1 + p{|+32*2| 2 + |*i^ 4 | 2 } 


= det 


( 53 ) 


□ 


B. Proof of Theorem 0 

For the single intermediate relay case, we use the equivalent 
representation 

Y = HX + W ( 54 ) 

of the scheme i.e. 


yi 53 


fti 0 


CO 

.52 54 _ 


_ 52*2 hi 


x 2 X4 


Wl w 3 

h 2 v 2 p + w 2 h 2 v 2}3 + W4 


= [1 + p\hi\ 2 ) [l + p{\Ag 2 h 2 \ 2 + IKAl 2 } 

- p 2 (\hi\ 2 \Ag 2 g 2 \ 2 ) 

= 1 + p\hi\ 2 + A 2 ^p\h 2 g 2 \ 2 + p|*-i | 2 + p 2 |*i| 

This gives that 

Pout(r) 

= Pr^log[l + p\hi\ 2 +A 2 (p\h 2 g 2 \ 2 + p|*-i | 2 + p 2 |*-i | 4 )] 

< 2 r log(p) 

and then for 

*1 := |*i| 2 ; *2:=|*2| 2 5 2 : = l52| 2 , 


( 55 ) 



























































the probability of outage is given by 


where 


' {G A2 ,g 2 }trO 


e -l h i+ h 2+92] d h 1 dh 2 dg 2 


{hi, h 2 , g 2 } £ O 

<^> log[l + phi + A 2 (ph 2 g 2 + phi + p 2 h\)] < 2rlog(p) 

i.e., 

{hi, h 2 ,g 2 } £ O -£=> 

1 + phi + A~(ph 2 g 2 + phi + p 2 h^j < p 2 *. 

At this point, we can see that since 

hi,h 2 ,g 2 < p° 

(from [5]), we can approximate outage as 

{hi, h 2 , g 2 } £ O => 1 + phi + ph 2 g 2 + phi + p 2 h l < p 2 ' 

and can thus write 

P out = [_ _ e-^ 1+ ^ + ^dhidh 2 dg 2 . 

d {hi,h 2 ,g 2 }€iO 

It follows that the effective (whitened) channel H e ft has the 
same outage as does the original channel H, and this outage 
curve is met by the use of the approximately universal 2 x 2 
CDA code. For the general case of having n — 1 intermediate 
relays, the whitened equivalent channel described by G33 is 
exactly given by This channel’s outage curve was found 
in [5], shown here in Proposition^ and is exactly met by the 
corresponding approximately universal 2 (n — 1 ) x 2 (n — 1 ) 
CDA code. □ 


Appendix IV 

Distributed perfect space-time codes 
A. Distributed perfect space-time codes 

As shown in [9], the basic elements of a CDA space-time 
code are the number fields F, L, with L a finite, cyclic Galois 
extension of F of degree n. For a being the generator of the 
Galois group Gal(L/F), we let z be some symbol that satisfies 
the relations £z = za(£) V ( £ L and z n = 7 
for some ‘non-norm’ element 7 £ F* such that the smallest 
integer t for which 7 * is the relative norm N- L / ¥ (u ) of some 
element u in L*, is n. From there we construct a cyclic division 
algebra D = L ® zh ® ... ® z ra_ 1 L. A space-time code 
X can be associated to D by selecting the set of matrices 
corresponding to the left-regular representation of elements 
of a finite subset of D. For an arbitrary choice of integral 
basis {/3i}™Zo for a submodule in the ring of integers Ol 
of L over F, the elements of the signaling set are of the 
form li = X)”=o fijfij, G 1, £ L, f iti £ Op. As a result, 
the CDA code-matrix form is 

n —1 / \ n 

X = Y. T ° ( dia 9 (4 • G ) ) = E (58 > 

j— 0 ' u =1 

with / = [fj.o fj,i ■ ■ ■ fj,n- 1 ] an d the concatenated / now 
being the n 2 -length sequence of Of information elements 


leaving the transmitter node (after normalization). Each n 2 x n 
matrix A u (corresponding to relay u) is created by first letting 

A n ,i = diag[ (3i cr(/3;) cr 2 (/3i) ••• cr ” _1 (A) ] nx „ 

for i = 0,1, ■ ■ ■ , n — 1 and then recursively creating A u l = 
Y n ~ u A n 7 , u = 1, 2, ■ ■ • , n where 

^°(/3o) ••• <r n -\(3 0 ) 


G 


r 


cr°{Pn-l) ■■■ tJ n 1 (/3 n _l) 


0 

1 

0 •• 

0 •• 

■ 7 
• 0 


A u ,o 

0 

1 •• 

• 0 

A u — 






Au,n— 1 


(59) 


(60) 


X 


Mi 

fA n 


(61) 


The perfect code requirement that the lattice generator ma¬ 
trix G and the power sharing matrices F- 7 , j = 0,1,• • • ,n— 1 
in (1601 . be unitary matrices guarantees that = I, 

essential for the code’s information losslessness. 
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